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; Abstract 

In this paper, we study a scalar conservation law that models a highly re-entrant 
manufacturing system as encountered in semi-conductor production. As a generalization 
of [^, the velocity function possesses both the local and nonlocal character. We prove 
the existence and uniqueness of the weak solution to the Cauchy problem with initial 
and boundary data in . We also obtain the stability (continuous dependence) of both 
the solution and the out-flux with respect to the initial and boundary data. Finally, we 
prove the existence of an optimal control that minimizes, in the L^-sense with p G [1, oo), 
the difference between the actual out-flux and a forecast demand over a fixed time period. 
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ufacturing system. 



1 Introduction and main results 

In this paper, we study the scalar conservation law 

pt{t, x) + {p{t, x)A(x, W{t)))^, =0, t > 0, < x < 1, (1.1) 

where ^ 

W{t) = / p{t,x)dx. 
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We assume that the velocity function A > is continuous differentiable, i.e., A G (^^([0, 1] x 
[0, oo)), in the whole paper. For instance, we recall that the special case of 

was used in 23 1. 



This work is motivated by problems arising in the control of semiconductor manufacturing 
systems which are characterized by their highly re-entrant feature. This character is, in 
particular, described in terms of the velocity function A in the model: it is a function of 
the total mass W{t) (the integral of the density p). As a generalization of (in which 
A = X(W{t))), here we assume that the velocity A varies also with respect to the local position 
X, as can be naturally encountered in practice. These phenomena also appear in some biologic 



models (modeling the development of ovarian follicles, see [17, |T^]) and pedestrian flow models 
(see I, I, ID). 

In the manufacturing system, with a given initial data 

p{0,x) = po{x), 0<x<l, (1.2) 

the natural control input is the in-flux, which suggests the boundary condition 

p{t,0)X{0,W{t)) =u{t), t>0. (1.3) 

Motivated by applications, one natural control problem is related to the Demand Tracking 
Problem (DTP). The objective of DTP is to minimize the difference between the actual out- 
flux y{t) = p{t, 1)X{1, W{t)) and a given demand forecast yd{t) over a fixed time period. An 
alternative control problem is Backlog Problem (BP). The objective of BP is to minimize the 
difference between the number of the total products that have left the factory and the number 
of the total demanded products over a fixed time period. The backlog of a production system 
at a given time t is defined as 

/3(t)= f p{s,l)X{l,W{s))ds- f yd{s)ds. 
Jo Jo 

The backlog f3{t) can be negative or positive, with a positive backlog corresponding to over- 
production and a negative backlog corresponding to a shortage. 

Partial differential equation models for such manufacturing systems are motivated by the 
very high volume (number of parts manufactured per unit time) and the very large number 
of consecutive production steps. They are popular due to their superior analytic properties 
and the availability of efficient numerical tools for simulation. For more detailed discussions, 
see e.g. ^, 23|. In many aspects these models are quite similar to those of traffic flows 



[^] and pedestrian flows [0, 10|. 

The hyperbolic conservation laws and related control problems have been widely stud- 
ied for a long time. The fundamental problems include the existence, uniqueness, regularity 
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and continuous dependence of solutions, controllability, asymptotic stabilization, existence 
and uniqueness of optimal controls. For the well-posedness problems, we refer to the works 



24, p^] (and the references therein) in the content of weak solutions to systems (in- 



cluding scalar case) in conservation laws, and [26, in the content of classical solutions to 
general quasi-linear hyperbolic systems. For the controllability of linear hyperbolic systems, 
one can see the important survey [Q. The controllability of nonlinear hyperbolic equations 
(or systems) are studied in [|^, U, while the attainable set and asymp- 

totic stabilization of conservation laws can be found in ||l|, 0]. In particular, provides a 
comprehensive survey of controllability and stabilization in partial differential equations that 
also includes nonlinear conservation laws. 

We prove the existence, uniqueness and regularity of the weak solution to Cauchy prob- 



lem (|1.1| ), (1.2) and ( |1.3D with initial and boundary data in The main approach is 

the characteristic method. We point out here that in the previous paper [^], the authors 
obtained the well-posedness for (1 < p < oo) data. The L°° assumption in this paper 
is due to the fact that the velocity function A depends on the space variable x. Using the 
implicit expression of the solution in terms of the characteristics, we also prove the stability 
(continuous dependence) of both the solution and the out-flux with respect to the initial and 
boundary data. The stability property guarantees that a small perturbation to the initial and 
(or) boundary data produces also only a small perturbation to the solution and the out-flux. 

The optimal control problem that we study in this paper is related to the Demand Tracking 
Problem. This problem is motivated by [jl5| and originally inspired by [^]. The objective 
is to minimize the L^-norm with p £ [l,oo) of the difference between the actual out-flux 
y{t) = p{t, 1)A(1, W{t)) and a given demand forecast yd(t) over a fixed time period. With 
the help of the implicit expression of the weak solution and by compactness arguments, we 
prove the existence of solutions to this optimal control problem. 

The main difficulty of this paper comes from the nonlocal velocity in the model. A related 
manuscript which is also motivated in part by [Q, |2^, addressed well-posedness for 
systems of hyperbolic conservation laws with a nonlocal velocity in M". The authors studied 
the Cauchy problem in the whole space M" without considering any boundary conditions and 
they also gave a necessary condition for the possible optimal controls. However, the method 
of proof and even the definition of solutions are different from this paper. Another scalar 
conservation law with nonlocal velocity is to model sedimentation of particles in a dilute fluid 
suspension, see [^] for the well-posedness of the Cauchy problem. In this model, the nonlocal 
velocity is due to a convolution of the unknown function with a symmetric smoothing kernel. 
There are also some other one-dimensional models with nonlocal velocity, either in divergence 
form or not, which are related to the 3D Navier-Stokes equations or the Euler equations in 
the vorticity formulation. Nevertheless, the nonlocal character in these models comes from 



a singular integral of the unknown function (see |16| and the references therein, especially 

iTii). 
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The organization of this paper is as follows: First in Section 2 some basic notations and 
assumptions are given. Next in Section 3 we prove the existence and uniqueness of the weak 
solution to Cauchy problem ( |1.1D , ( |1.2| ) and ( |1.3| ) with the initial data po S L°°(0, 1) and 
boundary data u G L°°(0,T). Some remarks on the regularity of the weak solution to the 
Cauchy problem are also given in Section 3. In Section 4 we establish the stability of the 
weak solution and the out-flux with respect to the initial and boundary data. Then in Section 
5, we prove the existence of the solution to the optimal control problem of minimizing the 
L^-norm of the difference between the actual and any desired (forecast) out-flux. Finally in 
the appendix, we give two basic lemmas and the proofs of Lemmas |3.1[ - p.2| that are used in 
Section 3. 



2 Preliminaries 

First we introduce some notations which will be used in the whole paper: 

1/^(0, 1) := {/ G L°°(0, 1) : / is nonnegative almost everywhere}, 
L'^{0,T) := {/ G L°°(0, r): / is nonnegative almost everywhere}, 
||Po||l°= := ||/Oo||l°°(o,i) := esssup |po(a^)|, 

0<x<l 

||u||loo := ||w||loo(q -7^) := esssup |n(t)| 

0<a;<T 

and 

M := \\u\\l^o,t) + WpoWl^o,!), (2-1) 

A(M) := inf X(x,W) > 0, (2.2) 

(x,iy)e[o,i]x[o,Af] 

ll^llco := l|A||co([o,i]x[o,M]) := sup |A(x,l^)|, (2.3) 

{a;,iy)e[0,l]x[0,Af] 

llA^rllco := l|Ax|lco([o,i]x[o,Af]) := sup \\^{x,W)\, (2.4) 

{a;,iy)e[0,l]x[0,Af] 

W^wWco ■= \\^w\\cO{[0,l]x[0,M]) ■= sup \Xw{x,W)\. (2.5) 

(x,M/)e[0,l]x[0,A/] 

We also define the characteristic curve ^ = (,{s; t, x), which passes through the point {t, x), 
by the solution to the ordinary differential equation 

^ = A(e(s),M^(s)), m = x, (2.6) 

where is a continuous function. The existence and uniqueness of the solution to ( |2.6D is 
guaranteed by the assumption that A G C^([0,1] x [0, cx))) with |VF(s)| < M for all s. The 
characteristic curve is frequently used afterward and it is precisely illustrated in different 
situations. 
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3 Well-posedness of Cauchy Problem with L°° data 

First we recall, from [jl2|, Section 2.1], the usual definition of a weak solution to Cauchy 
problem (0), Q and (|L3|) . 

Definition 3.1. Let T > 0, po £ L~(0, 1) and u G L°°{0,T) be given. A weak solution of 
Cauchy problem ^) and is a function p € C°([0, T]; L^(0, l))nL°°((0, T) x (0, 1)) 

such that, for every r G [0,T] and every ip £ C^{[0,t] x [0, 1]) with 

99(r, x) = 0, Vx G [0, 1] and (^(t, 1) = 0, Vt G [0, r], 

one has 

1 

x){ipt{t, x) + A(x, x))dxdt 

+ / u{t)(p{t,0)dt + [ po{x)(f{0,x)dx = 0. 



JO 



JO JO 

Theorem 3.1. Let T > 0, po G (0, 1) and u G L5^(0, T) be given, then Cauchy problem 
([Til), (|L2|) and (|L3|) arfmi<s a unique weak solution p G C°([0, T]; L^(0, 1)) nL°°((0,r) x 
(0, 1)), which is nonnegative almost everywhere in [0, T] x [0, 1]. Moreover, the weak solution 
p even belongs to C°([0, T]; LP(0, 1)) for all p £ [l,oo). 



Proof: Our proof is partly inspired from [|l^. We first prove the existence of weak solution 
for small time: there exists a small 6 G (0, T] such that Cauchy problem ( |1.1D , ( |1.2D and ( |1.3D 
has a weak solution p G C°([0, 5]; (0, 1)) n L°°((0,(5) x (0,1)). The idea is first to prove 
that the total mass W{t) exists as a fixed point of a map W i-^ F{W), and then to construct 
a (unique) solution to the Cauchy problem. 
Let 

^S,M ■■= \W G C0([0,5]): \\W\\co[o,s] ■■= sup \W{t)\ < m}, (3.1) 

0<t<S ^ 

where the constant M is given by ( p.lD . 

For any smaU 5 > 0, we define a map F : Q,s,M (^"([O, (5]), W -F(VF), as 

F{W){t):= / n(a)da+ / Po{/3)dp, VtG[0,5], (3.2) 







where (see Fig [T| or Fig ^) represents the characteristic curve passing through (t, 1) which 
is defined by 

^ = X{^,{s),Wis)), ei(i) = l- (3.3) 

Here we remark that the formulation of F{W) is induced by solving the corresponding linear 
Cauchy problem ( |1 . ![ ) , (1.2) and ( |1.3| ) in which VF(-) G ri^^A/ is known. It is obvious that F 
maps into Qs M itself if 



< 5 < min 



5 



Now we prove that, if 5 is small enough, F is a contraction mapping on Q.s^m with respect 
to the norm. Let W,W £ ^s,M and for any fixed t G [0, T], we define by 



Then, we have for every t € [0, S] that 
\FiW)it)-F{Wm\ = 

< IIpoIIl°° 



< ill/'oIlL"" • (^l|Ax||co|l?l - 6llcO([0,t]) + ll-^W^IIco||W^ - W^llcO([0,5]) 
By the definitions of and , we obtain 

ll?i-ei|lco([o,i])= \U0)-Cii9)\ 
o<e<t 



sup 

o<e<t 



< *ll^a;||c«llll -6llcO([0,t]) +i||AH/||co||W^ " W^llcO([0,<5]), 



thus 



Therefore, 



11^1 



lllcO([0,t]) 



< 



*ll-^H/||co I1TT7 



x\\c° 



■ Il^-^llc0([0,5])- 



\Fiwm-Fiwm\ < . iF-P7ibo(p,.]), vt E [0,5]. 



1 - t\\x 



xWco 



Let 6 be such that 



< 6 < min 



^ll'^xllco ' 4||/)o||l°° ll-^vyllc" ' ll^llco 



then 



(3.4) 



(3.5) 



\\FiW) - F{W)\\co^[o,s]) <^\\W- VF|bo([o,^]). 

By means of the contraction mapping principle, there exists a unique fixed point W = F{W) 
in flsM: 



W{t) = F{W){t) = / u{a)da + 
Jo 

Moreover, W is Lipschitz continuous: 



i~J^x{ii{e),w{e))de 



PoW)d/3, Vie [0,5]. 



W{t) = W{0) + [ W'{s)ds 
Jo 



with 



w'it) = u{t)-xi^iit),wit))poii- xi^i{e),w{e))de), te[o,6] 
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and thus 

I|W^'IIl«>(o,5) < + ||A||co||po||l°°- (3.6) 

Now we define the characteristic curve ^2 which passes through the origin (see Fig |l| and 
FigD by 

^ = X{C,is),Wis)), 6(0) = 0. (3.7) 

And then for any fixed t € [0,S], we define the characteristic curves S,3,Ca (see Fig|l| and Fig 
^) which pass through (t, x) by 

^ = mis), Wis)), Ut)=x, forxG [0,6(t)]. (3.8) 

^ = X{Us),W{s)), Ut) = x, for [e2(t),l]. (3.9) 

From the uniqueness of the solution to the ordinary differential equation, we know that there 
exist a G [0, 6] and (3 G [0, 1] such that 

6(a) = and 6(0) = /3. (3.10) 

Now we define a function p by 



u{a) 



-Jlx.iisW,wmde^ 0<x<Ut),0<t<6, 



p{t,x):={mW{a)) ' - -s.w, - - , ^^^^^ 

[po{f])e-fo^^iMe)me))de^ < ^^^^^ <x,0<t<6, 

which is obviously nonnegative almost everywhere in (0,(5) x (0, 1). Using the following two 



lemmas, we can prove that p defined by (3.11) is the unique weak solution to the Cauchy 
problem (0), (Q and ([Ol) . 



Lemma 3.1. The function p defined by ( |j.ll ) is a weak solution to Cauchy problem (1.1), 



( |1.2| ) and ( |1.3| ). Moreover, the weak solution p even belongs to C^{[0,6]; LP{0,1)) for all 
p € [1, 00) and the following two estimates hold for all t £ [0, 6]: 

0<W{t) = \\p{t,-)\\L^^o,i)<M, (3.12) 
||p(t,-)||L-(o,i) <e^"^^"^° -maxlllpoIlL-,^!^}. (3.13) 

Lemma 3.2. The weak solution to Cauchy problem ( |1.1[ ), ( |1.2| ) and ( |1.3D is unique. 
We leave the proofs of Lemma p.l| and Lemma |3.2| in Appendix. 



Now we suppose that we have solved Cauchy problem (1.1), {^^) and ( |1.3| ) to the moment 
T G (0,r) with the weak solution p G C°([0, r]; LP(0, 1)) n L~((0,t) x (0,1)). Similar to 



Lemma |3.1| and Lemma we know that this weak solution is given by 
p{t,x) -- 



Po(/3)e-io ^Aim,wie))d9^ -f q < < x < 1,0 < t < r, 

^fi^e-/:A.(C3W,H^W)^^^^ else. 
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Moreover, the two uniform a priori estimates ( 3.12| ) and (|3.13| ) hold for ah t € [0,r]. Hence 
we can choose S € (0, T) independent of r such that ( ^l5| ) holds. Applying Lemma and 
Lemma |3.2| again, the weak solution p € C^{[0,t]; Lf{0, 1)), as well as estimates ( p.l2| ) and 
( 3.13| ), is extended to the time interval [r, t + 5] n [r, T]. Step by step, we finally have a unique 
global weak solution p G C°([0, T]; L?'(0, 1)) n L~((0,r) x (0, 1)). This finishes the proof of 



Theorem 3.1 



□ 



Remark 3.1. Let p be the weak solution in Theorem 3.1 and W £ C*'([0, T]) be the total 
mass function: W{t) = p(t, x)dx. Let ^i, ^2, '^3, ^4 and a, /3 be defined by ( |3.3D , (^), (|3.8D, 
(|3.9| ) and ( |3.10| ) , respectively. It follows from our proof of Theorem 3T that (see Fig |^ Fig 
I and Fig |) 



p{t,x) 



PQ{l3)e--fo^-(Me),wie))de^ 



u{a) 



{X{0,W{a)) 
and the following estimate holds: 



if < 6(0 < x < 1,0 < t < T, 
else. 



(3.14) 



ll/>(t,-)llL-(o,i) <e'^"^^"^° •max{||/5o||L-,|^}, Vt G [0,r]. 



Moreover, W{t) can be expressed as (see Fig |5| 

u{a)da + 
u{a)da 



W{t) = 
which implies again that 



0^ JO 



\{M) 
) 



(3.15) 



(3.16) 



o<w{t) = \\pit,-)\\mo,i)<M, vtG[o,r]. 

Finally, W is Lipschitz continuous: 

W{t) = W{0) + / W'{s)ds, 
Jo 

where (see Fig ^) 



(3.17) 



W'{t) = < 



n(t)-A(6(t),mt))po(l- / X{Ci{e),Wi9))de), 



^ ' ^ ^^A(0,W^(er\0))) ^2 V ; - - 



and 



. C ||w||^oo 7^||\ II 1 

1^ l|L°°(o,r) < + ||A||c-o ■ max|||/)o||Loo, e " -iicoj.< 



oo. 



(3.18) 
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Remark 3.2. (Hidden regularity.) From the definition of the weak solution, we can expect 
p € L°°((0,T) X (0,1)) = L°°(0, l;L°°(0,r)). In fact, under the assumptions of Theorem 
3.1, we have the hidden regularity that p € ^([0, 1];LP(0,T)) for all p G [l,oo) so that the 



function t ^ p{t,x) € LP(0,T) is well defined for every fixed x € [0,1]. The proof of the 
hidden regularity is quite similar to our proof of p G C^{[0,T]; L^{0,1)) by means of the 
implicit expressions (|1|) for p and (|1|) for W{t) (see also ( ^.181) when T is large). 

Remark 3.3. If po G ^^([0, 1]) and u G ^^([0, T]) are nonnegative and the compatibility 
condition is satisfied at the origin: 

u{0) 



Po(0) = 0, 



where VF(0) = pQ{x)dx, then Cauchy problem ( |1.1D , (1.2) and ( |1.3D admits a unique 
nonnegative solution p G C°{[0, T] x [0, 1]). If, furthermore, po G C^{[0, 1]) and u G C^([0, T]) 
are nonnegative and the compatibility conditions are satisfied at the origin: 

u{0) 

A(0,IF(0))-^°(°^ = °' 

^^^"^'^°'^^"|A(0,t!o)Tp°''^^°^^'^^^°^ + ^(«))^o(0) + A.(0, WmPoiO) ^ 0, 

where W{0) = po{x)dx and W'{0) = u{0) - po{l)X{l,W{0)), then Cauchy problem ([l1|), 
(|1.2|) and (|1.3|) admits a unique nonnegative classical solution p G C^([0, T] x [0, 1]). 



4 Stability with respect to the initial and boundary data 

In this section, we study the stability (or continuous dependence) of both the solution p 
itself and the out-flux y with respect to po and u. That is to say: if the initial and boundary 
data are slightly perturbed, are the solution p and the out-flux y also slightly perturbed? 

Let J) be the weak solution to the Cauchy problem with the perturbed initial and boundary 
conditions 

^Pt{t,x) + {p{t,x)X{x,W{t))y^ = 0, t > 0,0 < X < 1, 
p{0,x) =po{x), 0<x<l, (4.1) 

p(t,0)A(0,Ty(t)) = n(t), 0<t<r, 

where W{t) := J^'p{t,x)dx. We denote that y{t) := p(l, t)A(l, VF(t)). We also define the 
corresponding characteristics with respect to the perturbed solution p: (as (|3.3D ), ^2 (^s 
I), (^3, a) (as dJ) and ( |3loD ) and (^4,^) (as and (p^)), respectively 



First we have the following theorem on the stability of the weak solution p. 
Theorem 4.1. For any e > 0, p G [1, 00) and any K > {) such that 

IIpo||l°°(o,i) + II'"IIl°°(o,t) < K, ||poIIl°°(o,i) + II'"IIl°°(o,t) < (4.2) 
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there exists rj = ri{£,p,K) > small enough such that, if 



\Po - Po||lp(o,i) + \\u - u\\lp{q,T) < V, 



then 



Mtr)-p{t,-)\\LP{o,i)<e, yte[o,T]. 
Proof: Solving Cauchy problem (^]^), we know from ( |3.17| ) and ( [4. 21 ) that 

< Wit) < IIpoIIlmo,!) + II^^IIli{o,t) <k, ytG [o,r]. 



(4.3) 
(4.4) 

(4.5) 



Replacing M by in the definitions of A(M) and ||A||(^o, HA^^Hco, ||Avy||(70 (see Section 2), we 
introduce some new notations as X{K) and (still) ||A||(70, ||Aa;||(70, HAvyHco in this section. 

We first prove the stability of the weak solution for small time. Let 6 be chosen by (3.5). 
For any fixed t € [0,5], we suppose that ^2(t) < ^2(0 (the case that S,2{'t) > S,2i't) can be 
treated similarly). In order to estimate ||p(t, •) — /o(t, ■) IIlp(o,i) for P S [IjCo), we need to 
estimate Jq^^^^ \'p{t,x) — p{t,x)\^dx, f^^^^^ |p(^i2;) — p{t,x)\^dx and /^^(i) \p{t,x) — p{t, x)\Pdx, 
successively. 



For almost every x G [0, .^2(i)]5 we know from Remark 3.1 that (see Fig 



\p{t,x) -p{t,x)\ 

^(") /^A,(«3(e),W{e)) de 



u[a} 



X{0,W(a)) 

< l^(«)^^(«)l ^-r^A.(f.w.me))c 

~ X{0,W{a)) 
\u{a)\ 



+ 



X{0,W{a)) 
u{a) 



X{0,W{a)) A(0,W^(a)) 



+ 



'.CxAi3{e),wie))de 



X{0,W{a)) 

< C\u{a) - u{a)\ + C\u{a)\\W{a) - W{a)\ 



+ c\u{a)\ / x^{C3{0),w{e))de 



X^{U0),Wi9))d9 

Here and hereafter in this section, we denote by C various constants which do not depend on 
t, X, p,]) but may depend on p and K. 

For the given u G L°°(0,T), let {un}n=i C C'HlO,^]) be such that n„ ^ u in LP{0,T). 
And for the given A G C^i[0,l] x [0,oo)), let {vn}^=i C C\[0,1] x [0,K]) be such that 
Vn — > Xx in C'^([0, 1] X [0,A']). Using sequences {un}'^=i, {vn}^=i and (|3.6|) , we obtain for 
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almost every x G [0,S,2i't)] that 
\p{t,x) - p{t,x)\ 
< C\u{a) — u{a) \ + C\un(a) — u(a) \ + C|n,„(a) — u{a)\ 
+ C\un{a) - Un{a)\ + C\W{a) - W{a)\ + C\W{a) - W{a)\ 

+ c\u{a)\ f \vn{lm^w{e)) - \^{l,{e),w{e))\ de 



+ c\u{a)\ nvn{UO),w{e)) - K{uo),w{e))\d9 



+ C\u{a)\ _vn{U0),w{e))d6- vn{U0),wi6))de 

< C\u{a) — u{a) \ + C\un(a) — u(a) \ + C|n„(a) — u{a)\ 
+ C'l^^(a)lll'^n - ^x\\cO{[o,i]x[o,K]) + Cn\a- a\ + Cn\u{a)\\a - a\ 
+ C„|n(a)|||e3 - e3lbo([s,t]) + CnHa)\\\W - M^||co([o,5])- (4.6) 

Here and hereafter in this section, we denote by C„ various constants which do not depend 
on t, 3;,/9,p but may depend on p^K and n (the index of the corresponding sequences, e.g. 
{un}n=i^ {vn}n=i and SO on). 

A similar estimate as ( |3.4D gives us 



lies - mc^im) < /"^r^'T -W- H^||co([o,5]). (4.7) 



From the fact that 

ft 



Ut)= [_xMe),w{d))de = x = ^3{t)= [ x{U9),w{6))de 

J a J a 

and the definition of X{K), we get also that 

\a-a\<=^^l\{U0),W{e))d9 

''\x{Uo),w{e))-x{uo),w{d)))d9 



1 

< J— ( t\\X.Jc" ll?3 - 6lb"(M) + t\Mco \\W - W^llcO([0,5]) 



A(A') 



< 



*ll-^W/||co 



\\W-W\\co^[o,5])- (4.8) 



A(A')(l-i||A,||co) 

On the other hand, by (|3.16D and Holder inequality, we have for every t € [0, 6] that 

\W{t)-w{t)\ 

{u{a) — u{a))da + / pQ^a) da — / pQ{a) da 

-'o -'0 

< s'^\\u-uhp^o,T) + Wpo - pohno,!) + IIpoIIl- [ m,ie),wie)) - x{^,{e),wie))\ de 
1 _ _ 

< 6i \\u-u\\lp(o^t) + WPo - Po||lp(o,i) 







+ ^l|Po||L-(|lAx||co||ei -6llc«([0,t]) + \\Xw\\c4W -W\\cO(^io^s])), 
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where q satisfies ^ + ^ = 1. By the definitions of and Cii we still have (|3.4|), thus 

\\W-W\\co^[o,s])= sup \Wit)-Wit)\ 

te[o,5] 



and furthermore, from the choice ( |3.5D of 6, 

\\W - T^||co([o,<5]) < Clin - n||iP(o,T) + C'll^o - Po\\lp{o,i)- (4-9) 



Therefore, combining ( |4.6| ), ( |4.7| ), ( |4.8| ) and ( |4.9| ) all together and using Lemma A.l, we 
obtain easily that 

K2(t) 



\p{t, x) — p{t, x)\^dx 

< C\\Un - 'w||^p(o,T) + ^11"^™ ~ '^^\\cO{[0,l]x[0,K]) 

+ CJ 

For almost every x € ['?2(0> 1]> '^^ know from Remark 3.1 that (see Fig |^ 
\p{t,x)-p{t,x)\ = |po(^)e-/o^^(?4W,WW)de _^^(^)g-/*A.(€4(e),iyW)d^^ 

<C\-p^(^)-pom\ + C\po{P)-pom 



For the given po G ^°°(0, 1), we let {p[f}^=i C C^lO, 1]) be such that ^ po in ^^(0, 1). 
With the help of sequences {PojJ^Li) {'Vn]'^=ii we obtain for almost every x G [^2(0) 1] that 

|p(t,x) -p(t,2;)| 

< C\-pS) - + C\pm - PoiM + C|p^(/3) - po(/5)| + C|/>^(^) - />^(/3)| 

+ c7|po(/3)| / |^n(?4W, w^W) - K{UO).w{e))\ de 



Jo 

+ C\pom [ \vn{UO),Wi0)) - XAU0),W{9))\d9 
Jo 

+ c\pom [ \vn{Uo),w{e))-vn{UiO),wie))\d9 

Jo 



< cipo(/3) - Pom\ + c\pup) - pom + c\p"om - pom + c\pom\\vn - a.iico([o,i]x[o,ai) 

+ a|^-/3| + C„|po(/3)|||e4-e4||cO([o,t])+^^n|po(/3)|||W-^||co([o,5]). (4.11) 



Similar to (3.4), we have 



11^4 - C4|bo([o,t]) < il^jl'iy, • 11^ - WWcoim). (4.12) 
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and in particular, 



1/5 - /3| = 1^4(0) - ^4(0)1 < , -W- W^llcO(M)- 



(4.13) 



1 — i II Aj: 11(70 

Therefore, by the choice (U) of (5, together with estimates (|49|) , (|4.1lD , (|4.12| ), ( [4.13|) and 
Lemma |A.2| , we get immediately that 

|p(t, x) — p{t, x)\^dx 

hit) 

< C'IIpo - Po||^p(o,i) + C'll^n - -^^'■Ilc'0([0,l]x[0,ii']) 



+ Cn\\ 



LP(0,1)' 



(4.14) 



Now it is left to estimate /^^^(j) l/^C^i^;) — p(t,x)\Pdx (see Fig|8|). Obviously, we have 

\p{t,x) - p{t,x)\Pdx < / |p(t,x) -p(t,a;)|Pdx+ / |^(t, x) - /9(t, x)|P(ix. (4.15) 



Then similar to estimates ( d.lOl ) and ( |4.14D , we get easily from (|4.15| ) that 

r?2(t) _ 

\p{t, x) — p{t, x)\'^dx 



6W 



< C'IIPo Po||^p(o,i) + C'll^n ^^llip(o,T) + ^11""" '^2;||^0([o,i]x[0,_ft']) 

^ ll_ lip ^ ||_ lip 

+ '^n\\U ~ U\\ i^p(fj^j'-^ + '-^n IlPo ~ PO II LP(0,1) ' 



(4.16) 



Summarizing estimates ( 4. 10 ), ( 4. 14 ) and ( 4.16 ), finally we prove for any fixed t € [0,(^] 
with 5 satisfying (p. 5]) that 



||p(t,-)-p(t,-)||^P(0,i) 
< C'IIPo ~ Po|lip(o,l) + C'll^^n - ^IIlp(0,T) + ^W'^n - ^a; ||^0([o,l] x [0,i^]) 
+ Cn\\u — "Ull^pj-Q^j.-) + CnllpQ — Po||/,p(o,l)- 



(4.17) 



Now if we take n large enough and then rj in (|4.3| ) small enough, the right-hand side of ( [4.17 ) 
can be smaller than any given constant e > 0. 

In order to obtain the global stability ( [4. 41) from ( [4.17| ), it suffices, according to ( |3.5D , 
to prove uniform a priori estimates for \\p{t, ||/)(t, •)||l°°(o,i), \\'p{tr)\\L^{o,i)^ and 

\\'p{t, •)IIl°°(o,i)- fact the desired a priori estimates are available by Remark |3.1| (see esti- 
mates (3.15) and ( |3.17| )), hence we finally reach ( [4.4[ ). This finishes the proof of Theorem 



4.1 



□ 



We also have the stability on the out-flux y. 



13 



Theorem 4.2. For any e > 0, p € [1, oo) and any K such that ( [4.2| ) holds, there exists 



rj = ri{e,p,K)) > small enough such that, if (4.3) holds, then 



\\y-y\\LP(o,T) (4.18) 



Proof: Fhst, we prove \\y — y||LP{o,<5) < e for (5 smah. Let 6 be such that (|3.5| ) holds. For 
any fixed t G [0,5], without loss of generality, we suppose that ^2(i) < C2(^) i^^^ case that 



^2{t) > ^2(0 can be treated similarly). By Remark 3J., we have for almost every t € [0,5] 
that (see Fig ^) 

\y{t)-y{t)\ 
= \pit,mi,wit))-pit,i)xii,wm 

< C\p{t, 1) - p{t, 1)1 + C\pit, 1)\\\W - W^||cO([o,5]) 
= C|po(?i(0))e"i'o - po(ei(0))e-/o 

+ C\p{t,l)\\\W-W\\co^[o,s]) 

< ciMlim - Po(ei(o))| + c\p{t, i)\\\w - Ty||co([o,5]) 

+ cipo(ei(o))i / (A,(ei(^),w(^))-A,.(ei(e),w^(^)))(i0 



JO 

Then using sequences {pq}^^i, {vn}'^=i, we obtain 

m)-y{t)\ 

< c\MUo)) - Po{Uo))\ + c\pmm - poium + cipucm) - po(6(o))i 
+ c\pmm - po(6(o))i + c\p{t, i)\\\w - M/|ico([o,5]) 

+ c|po(ei(o))| / \vn{li{0),Wi9))-vn{Ciie),w{e))\de 



Jo 

+ c|po(ei(o))i / K(ei(^),W(0))-A,(ei(^),Tr(e))i^i0 

Jo 

+ c\poiCim\ [ \vn{Ci{0),wi9)) - x.,{Cii0),w{e))\de, 

Jo 

< c|po(!i(o)) - po(!i(o))| + c|p^(|i(o)) - po{lM)\ 

+ C|p^(ei(0)) - Po(6(0))j + C|po(6(0))||bn - A.|bo([o,l]x[0,i^]) 

+ a|po(6(0))||ki-6llc«([o,t])+Cn|po(ei(0))||F-H^|bo([o,5]), 



which results in 
r<5 



/ \m-y{tWdt 

Jo 

^ C'II/Jq — Po\\^LP{0,l) ^11^" ~ '^^\\cO{lO,l]x[0,K]) + ^"ll/'O ~ ^'o||^p(o,i) + Cn\\u — '"||^p(o,T) 



with the help of ( |3.4[ ) and ( [4.91) . Obviously, |[y — y||LP(o,<5) < ^ holds if we let n large enough 
and then 77 > in (|4.3|) small enough. 
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In a same way, we can prove for every r € [0, T — 5] by applying Theorem 



4.1 



that 



\\y - y\\LP{T,T+s) < e 

if 6 satisfies ( |3.5| ) (independent of r). Step by step, we finally prove ( [4.18 ). 



□ 



5 L^-optimal control for demand tracking problem 



Let T > and po S (0, 1) be given. According to Theorem 3T, for every u € L'^(0, T), 
Cauchy problem (|L1|) , ([L^) and ( |1.3D admits a unique weak solution L°°((0, T) x (0, 1)) Pip € 
C70([0,r],LP(0,l)) nCO([0, 1],LP(0,T)) for all p € [l,oo). 

For any fixed given demand € L°°(0, T), initial data po G (0, 1) and p S [1, oo), we 
define a cost functional on L'^{0, T) by 

Jp{u) := ||'u||^cx)(Q^2") + \\y ~ yd\\\p{td,T)i ^ (0) 

where y(t) := p{t, 1)A(1, W{t)) is the out-flux corresponding to the in-flux u G LJ^(0, T) and 
initial data po- This cost functional is motivated by |l5| , and the existence of the solution 
to this optimal control problem is obtained by the following theorem. 

Theorem 5.1. The infimum of the functional Jp in L'^{0,T) with p € [l,oo) is achieved, 
i.e., there exists u°° G LJ^(0,r) such that 

JJu°°) = inf JJu). 
«eL^(o,T) ^ 

Proof: Let {ti"}^^ C L5|°(0,T) be a minimizing sequence of the functional Jp, i.e. 

lim Jp(n") = inf Jn{u). 

Then we have 

II^"IIl-(o,t) + l|y"llLP(o,T) < c, yn e z+. (5.1) 

Here and hereafter in this section, we denote by C various constants which do not depend 
on n (the index of the sequences, e.g. {n"}^^, {y^}^^i and so on). The boundedness of 
ti" C L^{0,T) implies that there exists € L'^{0,T) and a subsequence of {u'^'^}'^-^ such 
that n"'-' u°° in (0, T). For simplicity, we still denote the subsequence as {li^jj^]^. And 
we let 

M ■■= ||Po||li(0,1) + ™^^{ sup ||^i"||Li{0,T), II^^°°IIli(0,T)} < OO- (5-2) 
Let p" be the weak solution to the Cauchy problem 

pj'(t, x) + (p"(t, x)X{x, VF"(t)))^ = 0, t > 0, < X < 1, 
p"(0,x) =po(x), 0<x<l, 
p"(t,0)A(0,VF"(i)) =^x"(t), 0<t<r, 
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where W"'{t) := p^{t, x)dx. Then for any fixed t E [0, T], we define ^" by 



m{s),W''{s)), CiW = 1 for any fixed tG[0,r], 



and define ^2 by 



^ = A(C2"(s),Vr-(s)), e2"(0)=0. 



Thanks to (B.ll), we know from and ( p^ ) that 

II^"IIh/i.-(o,t) < VnGZ+, (5.3) 

Then it follows from Arzela-Ascoli Theorem that there exists € C'^([0,r]) and a sub- 
sequence {M^'^'l^j^ such that W"^^ — > in C''([0,T]). Now we choose the corresponding 
subsequence and again, denote it as {u"},^^. Thus we have vP' u°° in L'^{0,T) 

and VF" ^ W°° in CO([0,T]). 

In view of (p.l6|) , (|5.2[ ) and (|5.3|) , there exists a small 5 > depending only on M and 
independent of n, such that 

nt A(cj'(e),W"(e))d0 
= / u''(a)da+ / Po{P)dP, yte[0,6]. (5.4) 







For any fixed t S [0, T], we define by 



KKi{s),W {s)), ^1 (t) = 1 for any fixed tG[0,r], 



and define by 

A(e~(.),vF°°(s)), er(o)=o 



Obviously, in C°([0, T]) implies that for any t € [0, <5], ^ ?r in Ci([0, t]) and 

^2 — ^2° in C^([0,r]). Thus by passing to the limit n — )> cx) in ( ^.41) , we obtain 

W {t)= / n°°(a)(ia+ / /Oo(/3)d/3, VtG[0,(5]. (5.5) 



Let /9°° be the weak solution to the Cauchy problem 

' pf{t,x) + (p°°(t,x)A(x,Ty°°(t)))^. = 0, t > 0,0 < a; < 1, 
p°°{d,x) = pq{x), 0<a;<l, 
_p°°(t,0)A(0,VF°°(t)) = u°°(t), o<t<r, 

where W°°{t) := /J p°°{t,x)dx. For any fixed t G [0,r], we define by 
and define by 



A(C(s),VF-(5)), C(t) = l 



^ = A(e2~(^),VF-(s)), e2°°(0) = 0. 
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Let 5 > he small enough so that 



u°°{a)da + 



PoiP)dP, Vte[0,5]. 



(5.6) 



We know from the proof of Theorem |3.1| that W = F{W) has a unique fixed point in fi^ 
(replacing M by M in (|]|)). This imphes from (U) and (U) that W°°{t) = W°°{t) on 
[0,<5], hence for any fixed t G [0,5], Cf^(s) = (.T (s) , (s) = (.Tis) on [0,t]. Moreover, with 
the help of (^, there exists 5o > independent of r € (0, T) such that if W°°{t) = W°°{t) 
then W^{t) = TF°°(t) on [t,t + 6o] n [r,T]. Hence we have W in C°{[0,T]) and 

furthermore — > uniformly and ^2 ^ ^2° ™ C'^([0,T]). 

Next we prove that y"(t) = p"(i, 1)A(1, Vf^"(t)) converges to y'^{t) = p°°{t, 1)A(1, W°°{t)) 
weakly in (0, T) for p G [1, 00). By is bounded in (0, T). Hence, it suffices 

to prove that for any g G C'^{[0,T]), 

''^{y^{t)-y^{t))g{t)dt = 0. (5.7) 



lim 



Case i. if{T) < 1. 

r 

/o 



{y^it)-y^it))git)dt 



T 







,;9o(er(o)))A(i,Ty"(t))e-/o^^«rW'W^n^))'^^ 

-po(C(0))A(l,T^°°(t))e-/o^4«rW.^°°W)'^^)5(t)di . 
Let {Polfcli C C^([0, 1]) be such that ^ po m ^^(0, 1), then by ( |3l^ ), ( pD and Lemma 



A.2| , we have 



(y"(t)-2/~(t))g(t)dt 



< 



+ 



+ 



T 



T 



Po(^r(0))-po(er(0))|A(l,T^"(t))e-/o"^(«i('')'^"(''»'^^b(t)|dt 
Po(Cr(0))-po(Cr(0))|A(l,t^°°(t))e-i'oA.{?rW,H/-W)dej^(^)|^^ 

pg(^r(0))A(l,W^"(t))e-/o^^(«rW'^"(^))'^^ 

- /Oo(er(0))A(l,W'°°(t))e-i'o^--(«i°''('^)'^°°(^))'^'' |5(t)|dt 

< C\\pI - Po||lp(o,i) + Cu\\W^ - W^~|bo([o,T]) + Cfc r ler(O) - er(0)||5(i)M* 

JO 



(A,(er(e),i^"(e))-A,(cr(^),t^°°(e)))de i5(t)it^i 



(5.8) 



where is a constant depending on Pq. Noting the fact that — )■ po in L*'(0, 1) and 
^ ~^ uniformly, by Lebesgue dominated convergence theorem, we let k 
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large enough first and then n large enough, so that the right hand side of ( |5.8| ) can be 
arbitrarily small. This concludes ( ^.7| ) for the case of ^2°(r) < 1. 

Case 2. if{T) = I, i.e., T = (^f')"^!)- For every r G [0,r), we have 



< 



[y^[t)-y^{t))g{t)dt+ J {y^{t)-y<^{t))g{t)dt 
{y^{t)-y°-{t))g{t)dt +C{T-t)\. 



(5.9) 



Since it is known from Case 1 that for every r G [0,T), |/J'(y"(t) — y°°{t))g{t)dt\ ^ as 
n^oo, one has for T = (CD^^l) by letting r ^ T in (jsj) . 
Case 3. if{T) > 1. Now we have 

rT 



[ {y^{t)-y^{t))g{t)dt 
Jo 



+ 



{y^{t)-y^{t))g{t)dt. 



Using the result in Case 2, we need only to prove that | j'^^ooyi(^^^{y^{t) — y°°{t))g{t)dt\ — )• 0. 
For any fixed a G [0, T], we define by 



d^ 
ds 

d4r 

ds 



xi^^is),w-{s)), er(«) = o, 



And we know from W in C°([0,r]) that ^ and (^f)~^ (Cf)"^ uniformly. 

Then we get from (|1|) that 

{y^it)-y'^it))git)dt 



(«2°°)-Hl) 

T 



(p^it, 1)A(1, W^{t)) - p-(i, 1)A(1, W^mgit)dt 



n"(a")A(l,VF"(^)) ^_ |t„ x.{i-{9),w-{9))de 
A(0,W^"(a")) 

A(0,iy°°(a°°)) 

where a" = (Cf )~^(0), a°° = (^i°)~^(0) are determined by 



git)dt 



1 



and thus 



w''{9))de = 1 - / A(er(^), VF°°(0))(i0 = 0, 



da 



1 = A(i,w^"(0) A^crw.w'" 



A(0,VF"(a")) 
da°° ^ A(l,P^°°(t)) A.Kf 
dt A(0,H^°°(a°°)) 



(5.10) 

(5.11) 
(5.12) 
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Let (r", r°°) be the value of (a", a°°) when t = T in ( 5.10 ), and let (ry", 7]°°) be the value 
of (a",a°°) when t = (^f in ( ^l0| ). Since VF" ^ 1^°°, ^ uniformly, then 

r" — > r°°, r/" 77°° and a" — )■ uniformly. Hence, by (5.11)-( 5.12| ) and using the facts 
that n" ^ n°° and ^ ^f, (^f ^ (Cf")"^ uniformly, we obtain 

(2/-(t)-y-(t))<7(t)dt 



<C|77"-ry°°| + C|r"-T°°| + 



^"(«M(er)-'(i))-^°°(«M(er)"^(i)) 



da 



0, 



as n 00. 



This concludes the proof of ( [5. 7] ) for the case ^2^{T) > 1. 
As a result, 



Jp{u°°) — ||^°°|li=o(o,T) + 11?/°° ^rfllip{o,T) 

< liminf \\u"'\\Ioo(q t) + liminf - ydllipfoT) 

< liminf Jpfu") = lim JJu^) = inf J„(tt). 

This shows Uoo is a minimizer of Jp{u) in L!!p(0, T), and it proves also that — )• 7/°° strongly 
inLP(0,r). □ 



A Appendix 

A.l Basic Lemmas 

The following two lemmas are used to prove the existence of weak solution to Cauchy 
problem ( |1 . l]) , ( |1.2D and ( |1.3| ), when changing variables in certain integrals (see Section A. 2). 
We recall some assumptions that are given: A > 0, A € C^([0, 1] x [0,oo)), W € ^s.m (see 



(3.1) for definition). 



Lemma A.l. Let ,^3 he the characteristic (see ( |3.8D for definition) which passes through the 
point {t,x) and intersects the t-axis at the point (a,0). Then we have 



^ = ^ ^- ft\Aiizie),w{e))de 

dx X{Q,W{a)) 



(A.l) 
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Lemma A. 2. Let ^4 he the characteristic (see (|3.9| ) for definition) which passes through the 
point {t,x) and intersects the x-axis at the point (0,/3). Then we have 



dx 



(A.2) 



The proofs of Lemma A.l and Lemma A.2 are trivial and they can be found in [28|. 



The next lemma is useful to prove a uniqueness result (Section 2, Lemma 3.2) for Cauchy 
problem Q, Q and (|L3|) . 



Lemma A. 3. // p G C°([0, T]; Li(0, 1)) nL°°((0,T) x (0,1)) is a weak solution to Cauchy 
problem (1.2) and (|1.3|), then for every t E [0, T] and every ip € C-'^([0,t] x [0, 1]) such 

that 

(^(t,1)=0, VTG[0,t], 

one /las 

p{T,x){ipr{T,x) + X{x,W{T))(px{T,x))dxdT + / n(r)(^(T, 0)(iT 
Jo Jo 

1 .1 

p{t,x)ip{t,x)dx + / /3o(a^)'/3(0, x)(ix = 0. 





The proof of Lemma A. 3 is the same as the proof of Lemma 2.2 in |ll5| , Section 2.2]. 



A.2 Proof of Lemma 13.1 



First we recall some notations M, A(M), ||A||(70, IIA^Hco, ||Avi/||co which are defined in Sec- 
tion 2. And the characteristic curves Cii '^25 '?3, ^4 are defined by (|3.3|), (^^) , ( ^^ , (|3.9D, 
respectively. 

By definition ( p. 11 ), it is easy to see that p £ L°°{{0,5) x (0,1)) and that estimates 
( |3.12| ) and ( |3.13p hold. Next we prove that the function defined by ( p. 11 ) belongs to 
C*'([0, 6]; L^{0, 1)) for all p G [1, 00), i.e., for every t,t £ [0, 6] with t >t, we need to prove 

\\pit, ■)- P{t,-)\\LP (0,1) ^0, as l^-tj^O 

for all p G [1, 00). In order to do that, we estimate Jq^^^^ \p{t, x) — p{t, x)\Pdx, /^^(r) \p(^^ ^) ~ 
p{t,x)\Pdx and J^^l^f^^ \p{tj^) ~ p{t, x)\Pdx, successively. 



For almost every x G [0,£,2{t)], by ( 3.11 ) and noting W G ^s^m (see ( |3.1| ) for definition). 
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we know (see Fig |T^) 

\p(t,x) - p{t,x)\ 
u{a) 



u[a) 



X{0,W{a)) X{0,W{a)) 
u{a) 



X{0,W{a)) 
u{a) 



X{0,W{a)) X{0,W{a)) 
e 



X{0,W{a)) 

< C\u{a) - u{a)\ + C\u{a)\\W{a) - W{a)\ 

+ c\u{a)\ [\,{ue),w{e))d9 - [\,{U0),w{9))de 

J a J a 

where ^3 denotes the characteristic curve passing through {t,x): 

as 

and it intersects t-axis at (5,0): ^3(0) = 0. Here and hereafter in Appendix, we denote by 
C various constants which do not depend on x,t,t. 

For the given u G L~(0, T), we let {un}n=i C C^{[0, T]) be such that n„ ^ u in LP(0, T). 
And for the given A G Ci([0, 1] x [0,oo)), we let {vn}'^=i C Ci([0, 1] x [0,M]) be such that 
Vn — ^ Xx in C°([0, 1] X [0,Af]). Using the sequences {un}^=i, {vn}'^=i and noting (^), we 
have for almost every x £ [0,^2(0] that 

\p{t,x) - p{t,x)\ 
< C\u{a) - u{a)\ + C\u{a)\\W{a) - W{a)\ 



+ c\u{a)\ iXx{UO),w{e))de- Xx{im,w{9))de 

J a J a 

< C|u"(5) -n(5) I +C|'u"(a) -n(a) I +C|u"(5) -n"(Q) I + Cju(a) I |VF(5) -VF(a) I 

+ c\u{a)\ f \vn{UO), w{e)) - xSM.w{e))\de 

J a 

+ c\u{a)\ f \vn{UO), w{e)) - Xx{uo),w{e))\de 

J a 

+ C\u{a)\ fvn{U0),W{e))d9- [\n{U0),W{9))de 

J a J a 

< C|u"(5) - n(5)| + C\u"'{a) - u{a)\ + C\u{a)\\\vn - Ax||c"([o,i]x[o,M]) + Cn\a - a\ 

+ Cn\u{a)\\t- t\ + Cn\u{a)\\a - a\ + Cn\u{a)\\\C3 - 6llc"([5,t])- (A.3) 

Here and hereafter, we denote by C„ various constants which do not depend on x, t, t but 
may depend on n (the index of the corresponding approximating sequences, e.g. {un}^=i, 
{vn}'^=i and so on). 
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By the definitions of .^3,^3 and (3^), we derive 



Ms) -Us) 



and furthermore 



xiUe),wie))de - / x{U0),wie))de 

J s 

<C|?-t|+5||A,||co||e3-6llco([5,i]), ysG[a,t], 
U3-ei\\c0i[5,t])<C\t-t\. 



(A.4) 



Meanwhile, from the fact that 



Ut)= xiUO),wie))de = x = ut)= / xiUd),wie))d9 



and definition ( |2.2[ ) of A(M), we get 
1 



a — a < — 



A(M) 
1 



x{^3{e),w{9))de 



{Xiue),wi9)) - xiuo),wie)))de + / xiU0),wi9))de 



X{M) 

<CU3-U\c0i[a,t])+C\t-t\ 

< C\t-t\. 

Therefore, we get easily from ( |A.3| ), ( |A.4] ) ( |A.5| ) and Lemma A.l, that 

2{t) 



(A.5) 



x) — p{t, x)\^dx 

+ Cn\t-t\P. 

For almost every x & [S,2{t),l], hy definition ( f3.ll ) of p, we have (see Fig 11) 



<C\\Un uW^Lp^Q^^^ + C\\Vn "^2: ||^o ([0,1] x [0,Af]) 



(A.6) 



\p{t,x) -p{t,x)\ 
PO0) e-/o^^-«4(9),H/(e))de _ ^^^^^ MMe),wie))de 



< \po{^) - po(/3)| e-/^^(«^W'^W)'^^ + \pom 



<c\poi/3)-pom\ + c\pom / x.{U0).w{e))d9- x.{ue),w{9))de 

Jo Jo 
where the characteristic curve ^4 passing through (t, x) is defined by 

d^4 



ds 



X{U{s),W{s)), Ut)=x 



and it intersects x-axis at (0, /3): ^4(0) = /3. 
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For the given pQ G L°°(0, 1), we let C Ci([0, 1]) be such that po in (0, 1). 

Using sequences {/9g}J^^ and {vn}'^=i, we obtain for almost every x € [S,2it), 1] that 

\p{t,x) - p{t,x)\ 

< c\pu'^) - pom + c\pU(3) - pom + c\p^0) - pm\ + cipomiK - a.iico([o,i]x[o,a/]) 

+ C\pom\ fvn{U0),W{d))d9- f Vn{U0).W{e))d9 
Jo Jo 

< c\p^0) - po0)\ + c\p^m - pom + c\pom\K - a.iico([o,i]x[o,m]) 

+ Cn0 - /3| + C'nlpoimt- t\ + C„|po(/3)|||e4 " U\\cOi[0,t]) (^-7) 



By the definitions of 1^4,^4, we have 

IIC4 - U\\cO{[0,t]) < C\t - t\, 
which is similar to ( |A.4| ). In particular, 

|^-/3| = 1^4(0) -C4(0)|<C|t'-t|. 
Therefore, we get easily from ( |A.7] ), ( |A.8| ), ( |A.9 ) and Lemma A. 2 that 



(A.8) 



(A.9) 



6(0 



\p{t, x) — pit, x)\^dx 



< C\\po - Po|lip(o,i) + C\\Vn - Ax||^0([o,i]x[0,Af]) + ^'"1* ~ 



(A.IO) 



Finally, we turn to estimate /^^{t) l'^^*' ^) ~ P^^^ x)|^(ix (see Fig ^). Obviously, we have 

\p{t,x)-p{t,x)\Pdx < / \p{t,x)-p{t,x)\Pdx+ \p{t,x)- p{t,x)\Pdx. {A.n) 

Ji2(t) Jx 



Then similar to estimates ( [A.6| ) and (A.IO), we get easily from (A. 11) that 



\p{t, x) — p{t, x)\Pdx 



< C\\Un — u\\%(^o,T) + ^WPo ~ Po|Ilp(o,1) + C\\Vn — |[^o ([o,l] x [0,A/]) + ^"1* ~ (A. 12) 

Summarizing estimates (A^), ( A.10|) and ( A.12| ), we find that for every p € [l,oo). 



\\p{t,-) -P(*r)llip(0,i) 

< C\\Un — u\\'lp(^o,T) + ^WPo ~ Po\\^LP{0,l) + ^W'^^n — -^x ||^o (-[q^;^] ^ [0,Af]) ^'"1* ~ (A. 13) 

Therefore by Lebesque dominated convergence theorem, letting n large enough and then \t — t\ 
small enough, the right hand side of ( |A.13| ) can be arbitrarily small. This proves that the 
function p defined by ( |3TT| ) belongs to C°{[0, d];LP{0, 1)) for ah p € [1, oo). 
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Finally, we prove that p defined by ( p. 11 ) is indeed a weak solution to Cauchy problem 



(O), (O) and (O). Let r G [0,(5]. For any if G Ci([0,r] x [0,1]) with (^(r,x) and 



ip{t, 1) = 0, we have 

A:= I I p{t,x){(ft{t,x) + X{x,W{t))(fx{t,x))dxdt 



Jo 

r f<2{t) 



Jo A(0,T4^(a)) 



By Lemma A.l and Lemma |A.2|, we obtain 



A= r f u{a){y,t{t,Ut)) + KUt),W{t))^^{t,Ut))dadt 
Jo Jo 

+ j J Po{mMt,ut)) + HUt),w{t))^,{t,umdpdt 

Jo dt Jq Jq dt 

Hence by changing the order of integral, we get 



^0 Ja dt \ Jq Jq J fir) Jo ' dt 

where 

ft 



f{t):=l- / X{U0),W{9))de 
Jo 

represents the coordinate that the characteristic curve intersects with j;-axis and f^^{(3) 
represents the time when the characteristic curve staring from (0, /3) arrives at the boundary 
X = 1. Consequently, for any /? G [/(r), 1], and ^4 are identical to each other since they 
pass through the same point (0, /3), so we get immediately that 

and finally that 

A= r u{a){ip{T,Ur)) - v{a,Ua)))da + /^^"Vo(/3)(v?(r, ^4^) - V9(0, e4(0)))d/3 
Jo Jo 

+ / Pom^ir\f3),urHm - ^io,uo)))d(3 

Jf{r) 

= - / u{a)^{a,f))da - [ pQifS^O, (3)df3 - [ pQ{f3)ip{0, f3)df3 

Jo Jo J fir) 

= -[ u{a)ip{a,0)da - [ pQ{p)ip{0, P)dp. 
Jo Jo 

This proves that p given by ( 3.14| ) is indeed a weak solution to Cauchy problem ( |1.1D , 



rg) and (|TJ). 
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A. 3 Proof of Lemma |3.2| 

Let us assume that p G C°([0, 5];L^{0, l))nL'^((0, 5) x (0, 1)) is a weak solution to Cauchy 
problem ([1.1|), p.2| ) and ([1.3|). Here 6 is such that (3.5) holds. Then by Lemma |A.3| , for any 
fixed t £ [0, 6] and G C\[0, t] x [0, 1]) with (^(r, 1) = for r € [0, t], 



t rl 



p{T,x){(pr{T,x) + X{x,W{T))(px{T,x))dxdT + / n(r)(^(T, 0)(iT 



p{t,x)ip{t,x)dx + / po(a^)'/5(0, = 0. (A. 14) 



Jo 



where W{t) := p{T,x)dx. 

Let Vo e C^(0, 1), then we choose the test function as the solution to the following 
backward linear Cauchy problem 

^/;^ + A(x,TF(r))V'x = 0, 0<r<t, 0<a;<l, 
'4^{t,x) = i^o^x), < X < 1, 

^^(r,l)=0, 0<r<t. 

For any fixed t G [0,5], we define the characteristic curves CiiC2''?3i'?4 by 

^ = xcus),w{s): 
^ = xc^,{s),w{s): 
^ = x{Us),w{s): 
^ = x{Us),w{s): 

It is easy to see that there exist a £ [0,t] and /3 S [0,1] such that 

^3(0) = and ^4(0)=^. 

In view of ( A.14| ), we compute 

'p{t,x)-ipo{x) dx = / n(r)^(r, 0) (ir + / po{x)^l^{0,x) dx 





= 1, 




?2(0) 


= 0, 






= X, 


for xG [0,^2(0] 




= X, 


for xG [?2(i),l] 



u{-a)ii){a,{)) da + / po(/5)V'(0, /3) d/5 
;i-/o A(fi(e),W(e))de 

u{a)ilJo{x) da + 



Poil3)Mx)d/3 (A.15) 



By the definitions of a, /3 and similar to ( [A.1[ )-( |A^ ), we have 
dx X{0,W{a)) dx 



25 



Thus, ( [A. 15 ) is rewritten as 

x)ipo{x) dx 



\{0,W{a)) Jl^i^t) 

Since G C'o(O)l) and t G [0,(5] are both arbitrary, we obtain in C'^([0, 5]; (0, 1)) PI 
L°°((0,(5) X (0,1)) that 



p{t,x) -- 
which hence gives 



- !lx.me)me))<ie ^ < ^ < ^^(,) < 1,0 < t < 5, 

A(0, VF(a)) 

po(^)e-i'o^-(«4W,i^W)'^^', <?2(t) < X < 1,0 < t < 5, 



(A.16) 



= / -p{t,x)dx 
Jo 

rt ^i-/o'A(|i(e),W{e))<ie _ _ _ 

u{oi)d[a+ j po{l3)d(3 = F{W){t), VtG[0,(5]. 

'0 Jo 



By (3^), we claim that W G a/ and then W = W since is the unique fixed point of the 
map W I-)- F{W) in CIs,m- Consequently, we have = (i = 1,2,3,4) and a = a,/3 = /3. 
Finally, by comparing and (|Al^ ), we obtain p = p. This gives us the uniqueness of 

the weak solution for small time. 
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